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Abstract 

We study the problem to provide an implicit triangular form for non-linear multi- 
input systems with respect to the flatness property. Furthermore, we suggest a 
constructive method, how to transform a given system into that special triangular 
shape, if possible. The well known Brunovsky form, which is applicable with regard 
to the input to state linearization problem, can be seen as special case of this implicit 
triangular form. A key tool in our investigation will be the construction of certain 
Cauchy characteristic vector fields that additionally annihilate certain distributions, 
and that will correspond to variables whose time-evolution can be derived without 
any integration. 

Key words: Differential Flatness, Differential geometry, Pfaffian systems, 
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1 Introduction 

The concept of flatness introduced in [6,7] has greatly influenced the control 
and systems theory community. The property of a system to be flat allows for 
an elegant solution for many feed-forward and/or feedback problems and is 
applicable for a big class of systems including the linear and the nonlinear as 
well as the lumped- and the distributed-parameter case. Within this paper we 
are interested in the system class of nonlinear multi-input systems described by 
ordinary differential equations. For this system class necessary and sufficient 
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conditions for flatness have been proposed in [8,9] based on a polynomial 
matrix approach. Furthermore, for nonlinear multi-input systems with special 
structure further results exist, see e.g. [10,11,12,13]. 

Triangular forms are of special interest for nonlinear systems in the context 
of input to state linearizability or flatness. Systems that are input to state 
linearizable by static feedback can be converted to Brunovsky normal form (a 
special case of the extended Goursat form, adapted to control systems), see [5]. 
It is well known, that systems that are flat but not input to state linearizable 
by static feedback can be transformed into Brunovsky normal form only after 
a dynamic system extension (dynamic compensator), see [4]. For systems that 
are 0-flat in [2] a nonlinear explicit triangular form has been proposed. We 
consider systems that are 1-flat, the flat output may depend on the state (0- 
flat) and the control (1-flat) but not on the derivatives of the control. Our 
presented form relies on an implicit triangular decomposition, and we propose 
a constructive scheme how to transform a 1-flat system into that special form, 
if possible (this gives rise to a sufficient condition for a control system to be 
1-flat). As in [5] where the extended Goursat form is discussed, we also make 
use of the Pfafhan system representation, where 1-forms (covector-fields) are 
used for the description of the (implicit) differential equations. Furthermore 
we will also make use of a filtration, which will lead to an implicit triangular 
representation in the case of flat systems. This is in contrast to the filtration 
which is based on derived flags used in the input to state linearization problem 
leading to an explicit representation like the Brunovsky form, see again [5]. The 
implicit triangular decomposition contains the Brunovsky form as a special 
case. A different approach based on differential forms, but associated with the 
tangent linear system can be found in [1], where the so-called infinitesimal 
Brunovsky form is considered. 

This contribution can be seen as further developing the ideas presented in 
[14,15], where a reduction and elimination procedure is considered to derive 
flat outputs, but not based on a Pfaflian representation and the adequate tools 
from exterior algebra. Preliminary results have been presented already in [16], 
where also an extended example, the VTOL can be found, which has been 
analyzed using different tools e.g. in [7,8]. 

Notation: Let X be an n x dimensional manifold equipped with local coor- 
dinates x a , a = l,...,n x . We denote the partial derivatives by d x <* and 
d x = d x i , . . . , d x a , . . . , d x n x . We will make use of the Einstein convention on 
sums, namely aib 1 = J27=i a $ when the index range is clear from the context. 
If we have x = (x 1 , . . . , x Hx ) as well as x % = (x hl , . . . , x t,n ^) an expression of 
the form T = J2k=iYl7=i a ik( x )^ l ' k re& ds according to Einstein's convention 
as T = a J i k (x)x l ' k where each a,i(x) is a matrix of appropriate size and the 
indices j and k correspond to the rows and the columns of Oj(x). In matrix 
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notation we have 

a 3 ik (x)x l,k = ai(x)x 1 + . . . + a nx (x)x 
where each x % is identified with a column vector. 



2 The triangular form 

Let us consider a nonlinear control system 

x = f(x,u) (1) 

with n x states and n u independent inputs. Roughly speaking the system 
(1) is flat ((/t + 1)— flat), if there exist n u differential independent functions 
y(x, u,u, . . . , u^), such that the state x and the control u can be parametrized 
by y and its successive time derivatives. The system is called 0-flat if y depends 
solely on x and 1-flat if we have y(x, u). For a rigorous definition of differential 
flatness, see [6,8]. 

The Brunovsky-form, consisting of n u Integrator chains is the most simple 
triangular structure that can be achieved for systems (1) that are input to 
state linearizable by static feedback, and hence 0-flat. A different (explicit) 
triangular form for 0-flat systems has been proposed in [2] which is more 
general than the Brunovsky form. To treat the case of 1-flat systems we will 
present an implicit triangular form for a special subclass of systems of the 
form (1), that is useful regarding property to be 1-flat. 

Definition 1 The implicit differential equations S eji — , i — 1, . . . ,m (affine 
in the derivative coordinates z) 

- . Ij'i 1,1, ai _ lIJi 
— e,l - a l,ai Z 

■■ (2) 

^ . ™,jm - i,ai _ im,jm 
— e,m ■ u i,ai * u 

with ji = 1, . . . dim(S e i ) and oti = 1, . . . , n z i (i = 1, . . . , m) are termed the 
implicit triangular form (ji and are indices corresponding to the rows and 
the columns of the matrices a] and vectors b l , respectively) , which possesses 
the following properties 

(a) z = (z 1 , . . . , z\ . . . , z m ), where each z % is of the form z l = (z 1 ' 1 , . . . , z hn ^) 
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(b) z 1 = y 1 and we have a partition of the form 

z 2 = (y 2 ,z 2 ),...,z m = (y m ,z m ) 

with n z i = n y i + and n%i = 

(c) n z i > , i — 1, . . . , m , n y j > , j — 2, . . . , m 

(d) the matrices a\ and (the vectors) b % meet 

aZ=a^(z\...,z\^), (3) 
6 iJ ' , =6 i » J ' , (z 1 ,...,z i ,i i+1 ) 

(e) dim(S ej j) = dim(z i+1 ) 7 and the Jacobian matrices [d S i+iE eti ] are regular 
for all i = 1, . . . , m. 

Example 2 A system in triangular form with m = 3 in matrix notation reads 
as 

E^-.aKz 1 ^ 2 )^ -b\z\z 2 ) 

~ e ' 2 : a\{z\ z 2 , z^z 1 + a\{z\z 2 , z 3 )z 2 - b 2 (z\z 2 , z 3 ) 
E e>3 :a 3 (z\z 2 ,z 3 ,z 4 )z*-b 3 (z\z 2 ,z 3 ,z 4 ) 

where a summation from i = 1, . . . , 3 takes place and 

z 1 =y\ z 2 = (y 2 ,z 2 ), z 3 = (y 3 ,z 3 ) 

such that y 2 and/or y 3 are possibly empty (but they need not, n y 2,n y 3 > ). 
We require dim(S eii ) = dim(5 i+1 ) ; and the Jacobian matrices [d^+iE e ^] are 
regular for all i = 1 , . . . , 3 such that z l+1 can be computed by means of the 
implicit function theorem from S eji . 

Lemma 3 y is a flat output for the system (2). 

To proof this Lemma, we consider the implicit equations S ej i = 0. Then z l (t) = 
can be chosen freely, and z 2 (t) can be computed, where we make use of 
the implicit function theorem. We continue with the equations S e 2 = 0. Two 
scenarios are possible: z 2 = z 2 , then it can be easily checked that z 3 (t) can 
be computed using the same argument as for z 2 (t), since z 1 ^) and z 2 (t) are 
already given. If z 2 = (y 2 ,z 2 ) then y 2 (t) can be chosen freely, since the rank 
criteria is met for z 3 (t), which again can be computed. By continuing these 
procedure, we end up with the equations S e m = from which z m+1 {t) can be 
computed since at this stage z l {t), . . . , z m (t) are already known. This clearly 
shows that y 1 , . . . , y rn is a flat output for (2). 

Proposition 4 (a sufficient condition) The control system (1) is 1-flat if we 
can find locally a diffeomorphism (x, u) = y?(V, . . . , z m , z m+1 ) with n z i ) + 
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n^m+i = n u + n x , such that it can be represented in the form (2) and Y%Li n y' = 



The fact that this proposition is sufficient for 1-flat systems, comes from the 
observation that the flat outputs are among the z coordinates, and therefore 
clearly a function of x and special case also 0-flat systems are included. 

Furthermore it should be noted that we consider a diffeomorphism (x, u) = 
^(z 1 , . . . , z m , z m+1 ) which implies that we do not increase the dimension of 
the system variables. 

The goal is now to provide a constructive algorithm that transforms a non- 
linear multi-input system (1), if possible, into the form (2). Before we will 
analyze this in detail, let us consider an example. 



2.1 A motivating example 



We consider a system with 3 state variables (x 1 ,^ 2 ,^ 3 ) and 2 control inputs 
(u 1 , u 2 ) of the form 



x 1 — u 1 , x 2 = u 2 , x 3 = sin ( — | (4) 

\ u J 

also analyzed in [8,14] using a different approach. Let us introduce the local 
coordinate transformation (x 1 , x 2 , x 3 , u 1 , u 2 ) -H- (y 1 , z 2 , y 2 , z 3 , z 4 ) with 

x 3 



x 1 


= z 2 z 3 


y 1 


x 2 


= z 3 + y 2 


z 2 


x 3 


= y l 


y 2 


u 1 


= e z z z 


z 3 


u 2 


= e z 


z 4 



x 2 - x 1 ^ 

X 1 ^ 

II 1 - 



(5) 



\n{u 2 ). 

Then the system (4) in the new coordinates can be represented as 

y 1 — sin (z 2 ) = 

-y 2 z 2 + k 2 z 3 = (6) 

y 2 +z 3 -e i4 = 



(by a suitable combination of the equations) which is an implicit system of 
differential equations and the flat outputs are obviously y 1 and y 2 . Following 
Lemma 3 this can be seen as follows: Assigning we can compute z 2 (t). 
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Then by choosing y 2 (t) we can compute z 3 (t) from the second equation since 
z 2 (t) is already given. And from the last equation z 4 (t) follows since z 3 (t) 
and y 2 {t) are already given. Therefore, it is easily verified that the system 
variables (z 2 ,z 3 ,z A ) are parametrized with respect to the flat outputs y l ,y 2 
and its derivatives. Using (5) one obtains the system parametrization in the 
coordinates (x, u) and the flat outputs in original coordinates follow as y 1 = 
x 3 , y 2 = x 2 — x 1 ^. 

Example 5 The system (6) possesses the structure (2) as in example 2 with 
m = 3 and z 1 = y 1 , z 2 = (y 2 , z 2 ) , z 3 = z 3 and the matrices 

a\ = 1 , b 1 = sin (z 2 ) 

a\ = , a 2 = [-z 2 ,z 3 ] , b 2 = 

a 3 = , a 3 = [1 , 0] , a 3 = 1 , b 3 = . 

The key question is now, how to derive the coordinate transformation (5) 
and how must the equations be combined, such that the form (6) can be 
obtained. These two questions will be answered at once by using a Pfaffian 
system representation. 



3 Pfaffian representation 

We will use some facts about exterior algebra and Pfaffian systems in the 
sequel. For detailed information see [3]. More information concerning control 
systems in a Pfaffian representation can be found e.g. in [5] and references 
therein. It should be noted that we do not base our considerations on the 
tangent linear system, as it is used for instance in [1]. 

3.1 Introduction to Exterior Algebra 

We denote do; the exterior derivative of the k-form u and by v\ u the contrac- 
tion (interior product) of u by the vector field v. The exterior product (wedge 
product) is denoted by A. 

A Pfaffian system J on a n x - dimensional manifold X with coordinates x a is 
identified with a codistribution on X and is locally spanned by 1-forms ui on 
X. 

The derived flag of the Pfaffian System I is the descending chain of Pfaffian 
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systems J(°) D D /( 2 ) D . . .with /(°) = 7 and 



J (fe+1) := [uj G J (fc+1) ,du; = Omod/ (fc) } . 



Cauchy characteristic vector fields v of I denoted by v G C(I) meet 



(a)uJ/ = 0, (6) uJ(dJ) c/. 



(7) 



The importance of Cauchy characteristic vector fields lies in the fact, that the 
Pfaffian system I can be written using n x — c coordinates (after a suitable coor- 
dinate transformation), where c denotes the number of all independent Cauchy 
characteristic vector fields. Due to the fact that the distribution formed by all 
Cauchy characteristic vector fields is involutive, one can choose coordinates 
on X adapted to the vector fields, ), because of the 

Frobenius theorem, see [3], such that C(I) = {d x }. Because of (a) in (7) the 
forms in / do not include differentials dx and because of (b) in (7) one can 
chose a basis for I such that the coordinates x do not enter the function part of 
the forms, i.e. I can be expressed by the n x — c coordinates x, for the proof see 
again [3]. The same holds true for any involutive subdistribution C(I) C C(I) 
with dim(C(/)) = c < c such that I can be expressed by the n x — c coordinates. 

The annihilator of a Pfaffian system J is a distribution denoted by I L := 



3.2 Properties of time-invariant implicit Systems in Pfaffian representation 

To handle time-invariant implicit dynamical systems denoted by S we consider 
a fibration of a manifold Z over a one dimensional manifold (equipped with 
the time coordinate t), such that S is spanned by 1-forms reading as 



on a manifold with dimension n z + 1. To u ! there correspond the implicit 
differential equations u l e = with u l e = m % a (z)z a — n % (z). 

Definition 6 The vertical annihilator of (8) denoted by V '(S) 1 - consist of ver- 
tical vector fields (tangential to the fibration) w = w a (z)d z a that annihilate 



Example 7 Let us consider the control system x = f(x,u) written as a Pfaf- 
fian system S = {u) a } with 




uf = m^^dz - n l (z)dt 



(8) 



S. 



u a = dx a - f a (x,u)dt. 



Then we have S 1 - = {d t + f a (x, u)d x a, d u } as well as V(S) 



_L 



{d u }. 
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If the Pfaffian system S as in (8) takes the following special form 

w' = m l a (z, z)dz a - n l (z, z)dt (9) 

then the variables z are termed non- derivative variables. It should be noted 
that d z e V(S')~ L holds. If we furthermore assume that the implicit differential 
equations corresponding to (9), i.e. u\ = with oj l e = m l a (z,z)z a — n l (z,z) 
fulfill that the Jacobian matrix is regular and i = 1, . . . ,n z then we call 

the system S parameterizable with respect to z since then I = g(z, z), which 
follows at once by the implicit function theorem. 

3.3 The implicit triangular form in Pfaffian representation 

Let us consider the implicit differential equations as in (2) written using dif- 
ferential forms with the same properties as described in definition 1 (a) - (e) 
where Hj corresponds to the Pfaffian representation of 5 e> j 

E^.al'^dz 1 ' 01 -b^dt 

E 2 : alfldz 1 '^ + a 2 2 %dz 2 ' a * - b 2 ^dt 

'■■ (10) 

E m :a^ a ] m dz i ' a ' -b m ' lm dt 

Let us denote byLl] S^o the system (10) and by Sd,k — {Si, . . . , E m ^k}- 
Proposition 8 The system (10) enjoys the following properties 

(a) V(Sd,k) ± = {d sm +i-k} are involutive distributions and dgm+i-k G C(Sd,k+i) 
for k = 0, . . . , m — 1. 

(b) Each subsystem E^ is parameterizable with respect to the non-derivative 
variable z h+1 , i.e. 

z k+1 =g k+1 (z\...,z k ,z\...,z k ) 

for k = 1, . . . , m. 

(c) If in Ek, z k = (y k ,z h ) such that n z k > n z k, i.e. variables y k are present, 
then, d y k E C(S d , m - k +i)- 

The proof of this proposition is straightforward and follows from the structure 
of (10) together with the special structure of the ct^ 1 and 6 Mi according to 
(3) as in definition 1. 

1 The subscript d will always refer to a representation based on the desired trian- 
gular decomposition (10). 
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Corollary 9 The implicit triangular decomposition (10) gives rise to the de- 
composition of Sdfi into a nested sequence 



. . . c Sa,2 c Sd,i c Sd,o 



as well as to splittings of the form S^i = <Sd,i+i © <Sd,i+i,c where all the Sd,i+i, c 
are parameterizable, with respect to the corresponding non- derivative variables 
z. 

Example 10 (Example 2 cont.) Following the notations in proposition 8 we 
have S d ,o = {Si,S 2 ,S 3 } ; S d ,i = {Si,S 2 } and S d ,2 = {Si} since m = 3. From 
(a) we observe that V(Sd,o) J ~ = {d^} C C(Sd,i) which is obvious since z A 
is a non- derivative variable which only appears in S 3 . The same holds true 
regarding z 3 where now V(Sd t i) ± = {dp} C C{Sd,2)- Proposition 8 (c) means 
for instance that if z 2 = (y 2 , z 2 ) then d y 2 G C{Sd,2) since in Si only z 2 appears. 



4 A constructive algorithm 

The goal is now to develop a constructive scheme that subsequently creates 
this nested sequence based on a given control system of the form S = {cOq} 



The starting point of the scheme is the explicit system So but since linear 
combinations of the Uq lead to implicit equations in general we demonstrate 
the constructive method with the system Sk = {^1} (here the index k refers 
to the k — th iteration of the reduction process) with 



with i — 1, . . . , n e and n s > n e , where we denote by x all the system variables. 
(11) is a special case of (12), i.e. x = (x, u) in S . The following steps need to 
be performed 

(a) Computation of V(Sk) L , since these elements correspond to non- derivative 
variables. Choosing of an involutive T k C V(Sk)' L corresponds to a selec- 
tion of non- derivative variables called Wk- (This correspondence becomes 
obvious in an adapted coordinate chart to be constructed by means of 
the flow box theorem.) 

(b) Construction of a splitting Sk = SWi © «Sfc+i,c such that Tk C C(Sk+i), 
since this guarantees that Sk+i is independent of Wk- 

(c) Check, if Sfc + i )C is parameterizable with respect to the Wk, which is pos- 
sible only if dim(<Sfc) = dim(S , / t + i) + dim(J r fc ) holds. 

The whole procedure will then be continued with S^+i- 



uj« = dx a - f a (x,u)dt. 



(11) 



co k = m i a (x)dx a - n i (x)dt 



(12) 
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4-1 The k-th step of the system decomposition 

The constructive scheme rests on the following proposition 

Proposition 11 Let us consider the system S k = {u l k } with uu k as in (12). If 
we find an involutive distribution Tk withT k £ V(S k ) ± and a sub-codistribution 
Sk+i C S k such that T k G C{S k+1 ) is met, then we obtain a splitting S k = 
Sk+i © S k+1>c with 

Sk+i ■ = al{w k )dw k a - tf(w k )dt 

Sk+i,c ■ ^l+i,c = a ] ac {w k ,w k )dw k a - b> c (w k ,w k )dt 

in adapted coordinates (w k ,w k ) by using a diffeomorphism x = (p k (w k ,w k ) 
with n s = n Wk + n Ak . 

To prove this proposition we consider the case dim(J r fc ) = 1 first. A distribution 
generated by a single vector field v k is involutive by construction, and one 
performs a coordinate change based on the flow of v k which is denoted by <fiu> k 
with the flow parameter w k . If we change coordinates according to 

x = <pu, k ° ip(w k ) = <p k (w k , w k ) , n s = n Wk + , = 1 

where if) is chosen such that ip k is a local diffeomorphism, we obviously obtain 
in new coordinates v k — > d^ k and it can be checked, that d& k C V(S , / t)- L as 
well as dw k C C(S k +i) is met, therefore no dw k can appear and a basis of S k+ i 
must exist which is independent of the w k coordinates, since dw k C C(S k+ i). 

Remark 12 The flow of 4>w k is a ma P % — where x, corresponds to the 

initial conditions. Since the flow parameter w k is part of the new coordinates, 
and we require a diffeomorphism, the map x = ip(w k ) is at need, which cuts 
the flow (f>ti k transversal. 

To generalize this to dim(J-fc) > 1 is straightforward, since then the coordinate 
transformation is based on the composition of the flows generated by the fields 
of T k (where the pairwise Lie- brackets of these fields have to vanish). 

Based on these considerations we state the following corollary which addi- 
tionally includes the parametrization criteria, such that proposition (11) is 
connected with the triangular form (2), respectively (10). 

Corollary 13 The system S (11) can be transformed into the form (10) if 
we find a nested sequence of codistributions 

...cS 2 cS 1 cS 
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as well as involutive distributions Ti with dim(J-i) = l r that meet T\ G V(S'i)- L 
as well as T\ G C(Si + \) for I > such that the systems Si+i jC according to 
Si = Si+i © <Si+i, c are parameterizable. Then also dim(Si) = dim(5j + i) + l r 
holds. 

This sequence ends when we have a decomposition of the form Sk* = Sk*+\ © 
Sk*+i,c with Sk*+i the empty system, which means that Sk* is a parameter- 
izable system. This iterative scheme has therefore to be continued until a 
parameterizable system is obtained. It should be noted that in practice the 
effective computation of Ti and Si+i such that additionally parametrization is 
guaranteed for all elements of the sequence is a difficult task. We will comment 
on computational issues in the forthcoming and demonstrate on an example 
a possible strategy. 



4-2 The connection with the derived flag 

In this short paragraph we want to discuss how the derived flag, see [3,5] is 
connected to our filtration, as in corollary 13. 

Let us introduce an adapted basis for Sk, dim(S'fe) = h (with respect to the 
derived flag) which is Sk = {f^, . . . , £}? k , . . . , The first derived flag 
of Sk, denoted by S^\ meets Cl{ G S^ C Sk , j — 1, . . . , dim(S^) such that 

dti{ = of A ti[ + ft A Ql (14) 

holds for suitable 1-forms o>\,j3 J r . 

If {S^\dt} is integrable (Frobenius theorem, [3]), we have furthermore that 

dn{ = 7/ A ti[ + fP A dt (15) 
is met, for suitable 1-forms 7/,^. 

Proposition 14 Let us consider any sequence of Pfaffian systems . . . C S2 C 
Si C So with S as in (11) with elements Sk together with their first derived 
systems s£\ Then for all v G V(Sk) ± we have 

(a) vldS^ C is met if {Sj^dt} is integrable. 

(b) v\{dSP)cS k 

The proof of the first claim (a) follows by evaluating fjdf^ using (14) and 
(15) for v G V( y Sk) J ~ and the second (b) can be shown by using (14) in a 
straightforward manner. 
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Systems that are input to state linearizable by static feedback meet {Sj*\ dt} 
is integrable for every k, see [5]. 

Corollary 15 From Proposition 14 it follows that for systems S that are 
input to state linearizable by static feedback the sequence of the derived flags, 
corresponds to the sequence as in Corollary 13 and Tk corresponds to V(S k ) ± 
which is integrable by construction. 

Remark 16 The interesting case are of course examples that are not input 
to state linearizable by static feedback, i.e. {Sjp,dt} are not integrable, since 
then a different filtration has to be considered that may lead to an implicit 
triangular form. 

4-3 A constructive method to derive T k and S k+ i 

If one is able to construct the sequence as in Corollary 13, then one eventu- 
ally end up with the form (10) (by relabeling the coordinates) where in each 
step the involutive distribution J 7 / has to be integrated, in order to derive 
the coordinate transformation. Thus, in principle a constructive method that 
generates the implicit triangular decomposition is stated. If then the rank and 
dimension condition as in corollary 13 hold the system is 0-flat/l-flat, but it 
should be stressed, that this method is only sufficient for flatness, and a failure 
does not in general proof that a system is not flat. 

However, the construction of S k+ i C Sk such that an involutive distribution 
J-fc can be found that meets T k G V(Sk) J ~ as well as Tk G C(S k+ i) is a difficult 
task and lead in general to partial differential equations. Furthermore, since 
the choice T k G V(Sk) ± as well as of S k+ i with T k G C(S k+ i) is not unique in 
general (may lead to an 'dead end' of the algorithm) it might be necessary to 
iterate the construction of Tk and Sk+i (see section 5.2) - it should be noted 
that based on a simple necessary condition candidates for T k and Sk+i are 
singled out as shown next. 

We have to construct JF k G V(S k ) ± and S k +i C Sk such that T k \dS k+ i C S k+ i 
is met. Then also the necessary condition 

T k \dS k+1 C S k (16) 

holds, since S k +i C Sk- For V(Sk) ± = {vi} and S k = {w 3 '} with r = dim(S , fe) 
we derive the purely algebraic conditions (necessary conditions) 

c^Jd(aX) A (a; 1 A, . . . , Aco r ) = (17) 

where c 1 and aj depend on all the system variables, i.e. c l {x) and aj(x) when 
considering the system (12). 
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It should be noted that due to the requirement dim(S'fc) = dim(Sfc + i)+dim(J-fe) 
one has to find dim(iSfc) — dim(J r fc) independent solutions for S^+i- From 
(b) in Proposition 14 we see that fulfills this necessary condition indepen- 
dently of J-fc- If we furthermore assume that C Sk+i then the construction 
of Sk+i and Tk can be simplified further as demonstrated in the next section 
in great detail. The strategy is now to solve the necessary conditions (16) or 
which is the same (17) and to generate solutions for which then finally the 
criteria T^\dS^ + \ C Sk+i has to be checked as well as the parametrization as 
in corollary 13 . 



5 Examples 



We now present two examples, in the first one we show how one algorithmically 
can compute the nested sequence of codistributions using the necessary con- 
dition (17) and the second example demonstrates a case where the algorithm 
stops in a dead end, and another iteration is at need. 



5.1 The motivating example revisited 

Let us write the equations (4) as a Pfaffian system of the form So = {o;q, uj 2 , cOq} 
Uq = dx 1 — u l dt 

u 2 = dx 2 -u 2 dt (18) 
ujI = dx 3 - sin dt 

then we obtain the following proposition regarding the first reduction step 

Proposition 17 Given the system S as in (18) we derive a splitting of the 
form S = Si © Si >c as well as v that meets v e V(So) 1 " and v e C (Si) with 

vq = v}d u i + u 2 d u 2 

and Si = {u{ = Uq, uj\ = u 2 u)q — v}uJq) 

w\ = da; 3 - sin (^j dt , u\ = vfdx 1 - v}dx 2 (19) 

as well as a complement Si jC = {ujf c = uj 2 = dx 2 — u 2 dt}. 
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The proof of this proposition follows from the observation that V(Sq) ± = 
{d u i,d u z} and that voJdSi C Si as desired. We will now show how one can 
derive Si and v . 

Calculation 18 The first derived system Sq 1 ^ is given by the single form 

$ = cosf^rr) [ ^;dx 2 - dx 1 } + u 2 (dx 3 - sin(^r)df) 
u \u z J ir 

and a basis for S can be alternatively given as S = {uq, u 2 , 

To construct Si we assume that Sq^ C Si and consider the relation ( according 
to (17)) 

(c\d u i + c\d u 2)\d(a\ul + a^ 2 + a\$) AwjAw 2 A$ = (20) 
where c\ and a] are functions of all system variables, that have to be computed. 

2 

From (20) we are left with the equation c\a\ + c\a\ = or a\ = — a\^. This 

2 

means that the forms $ and u 2 — ^jUq fulfill the necessary conditions for the 

vector field v = c\d u i +cfd u 2. To determine c\ and c\ we consider the criteria 
Vi\(dSi) C Si and we derive the relation 

((cfoi + cjd u 2)\ (d$)) A K 2 - l^ 1 ) A $ = 0. (21) 

For the solution of (21) of the form c\ = c 2 ^ we have that $ and cJq — 
^tOJq clearly correspond to Si as in (19) as can be checked easily (by linear 
combinations) and that 

(v \ du 2 ) A u\ A u\ = (22) 
is fulfilled, such that u J(dSi) C Si is met, because of (21,22). 
The flow of vq, is given as (x, u) = 4>w{x, u) with x % = x % for i — 1, 2, 3 and 



u 2 = e™u 5 

with flow parameter w. We set x % = w l for i — 1, 2, 3, u 4 — w 4 and u 5 = 1 
(this corresponds the choice of tp, i.e. (x,u = ip(w)) to obtain locally the dif- 
feomorphism (x 1 , x 2 , x 3 , u 1 , u 2 ) = <£o{w l , w 2 , w 3 , w 4 , w) and we obtain a basis 
for S'i as 

uj] : dw 3 — sin (w 4 ) dt 

1 V ' (23) 

uj\ : dw 1 — w 4 dw 2 
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and for complement Si iC 



cul c : dw 2 - e*dt, 

and it can be checked easily that dim(J r ) = 1, dim (Si) + 1 = dim (So) and 
that the Jacobian d w (w2 — e w ) has maximal rank. 

Remark 19 We want to point out again, that v is a Cauchy characteristic 
vector field for S±, i.e. Vq G C(Si) and this guarantees that there is a basis 
for Si which does not depend on the coordinate w, since in new coordinates 
d w G C(Si) is met. 

Then we continue our considerations with Si and the following proposition 
states the second reduction step 

Proposition 20 Given the system Si as in (23) we derive a splitting of the 
form Si = S2 © S2, c as well as v± that meets v\ G V(Si) ± and v± G C(S2) with 



vi = w 4 d w i + d w 2 (24) 

and 

S 2 = = dw 3 - sin (w 4 ) dt} , (25) 
and S 2iC = {oj\ c : dw 1 — w^dw 2 } 

The proof follows again from the fact that V(Si)- 1 = {w 4 d w i + d w 2,d w 4.} and 
ViJdS 2 C S 2 . The construction of S 2 can be performed in the same manner as 
above. (Observe however that Sj 1 ^ is empty, but from 

(cl(w A d w i + d w 2) + c\d w i)\d{a\uj\ + a%uf) Aw}AwJ = 

that result follows at once.). 

Based on the flow of v\ we derive the map 



1 - 4 1 1 
w =qq + q 

w 2 = q + q 2 

w 3 = q 3 

4 4 
w = q 



with flow parameter q and we make the choice q 1 = (this corresponds again 
the choice of if)) leading to the transformation w = tpi(q, q). With 

(y 1 = q\ y 2 = q 2 , z 2 = q\ z 3 = q, z A = w) (26) 
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it is easily seen that the composition of (x, u) = <po(w, w) and w = <pi(q, q) to- 
gether with (26) gives the desired transformation (5). Furthermore the nested 
sequence of systems S 2 C Si C Sq leads at once to the desired normal-form 

^io = dy 1 ~ sin (^ 2 ) dt 

u d — zdz — z ay 
uj = dz 3 + dy 2 - e^dt. 

as in (6). The flow parameters w and q correspond to the non- derivative vari- 
ables z 4 and z 3 , respectively. 

5.2 A further example 

Let us consider the system So = {uq,Uq, Uq, u 4 } also treated in [4] in a different 
context 

= dx 1 — (x 2 + x 3 u 2 )dt 
u 2 = dx 2 - (x 3 + x l u 2 )dt 
uj 3 = dx 3 - (u 1 +x 2 u 2 )dt 
u 4 = dx 4 — u 2 dt 

where we again have V(So) 1 " = {d u i,d u 2}. 

The triangular form is based on the decompositions So = Si © Si )C with 

51 = {u\ = u^,u 2 = u 2 , uf = u 4 } 

uj\ = dx l - (x 2 + x 3 u 2 )dt 
uf = dx 2 - (x 3 + x 1 u 2 )dt 
u 3 = dx 4 - u 2 dt 

and Si :C = {<jj\ c = dx 3 — (u 1 + x 2 u 2 )dt] as well as on Si — S 2 © S 2jC with 

5 2 = {to\ — oj\ — u 2 uj, u 2 = uf} with 

u\ = dx l - u 2 dx 2 - (x 2 - x\u 2 ) 2 )dt. 
u 2 = dx 4 — u 2 dt 

and S 2>c = {oj 2c = dx 2 — (x 3 + x 1 u 2 )dt}. The distributions ~Fq = d u i G 
V(S'o) ± , Ti = d' x z E V(Si) 1 - and T 2 = d x 2+u 2 d x i e V(S , 2 ) ± were used and the 
flat outputs y 1 = x l —u 2 x 2 and y 2 = x 4 follow at once by applying a coordinate 
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transformation based on the flow of Ti and regarding JF$ and T\ no coordinate 
transformation is at need, since u 1 and x 3 are already non- derivative variables. 

Remark 21 Also in this example it holds true that Sq 1 ^ C Si and C S 2 
which enables one to construct the solutions based on the necessary condition 
(17) very easily. 

However, a different possible solution for Si © Si iC can be based on choosing 
the distribution Fq = {d u i,d u 2} together with Si = Sq 1 ^ = {uj\,uj 2 } 

uj\ = dx 1 - x 3 dx 4 - x 2 dt 
uj\ = dx 2 — x 1 dx i — x 3 dt 

and Si )C = {dx 3 — (u 1 + x 2 u 2 )dt , dx 4 — u 2 dt}, where obviously Si jC is param- 
eterizable with respect to u 1 and u 2 . This choice for T§ and Si however leads 
to an 'dead end' since for Si the necessary condition (17) does not lead to a 
splitting Si = S 2 © S 2iC . Therefore this choice of splitting has to be dismissed 
and one has to start again with S (iteration). 



6 Discussion 

We have characterized a suitable normal form for 1-flat systems, which is in 
implicit triangular shape, see (10) that possesses the properties as in proposi- 
tion 8 based on exterior algebra. Furthermore, we have discussed a construc- 
tive calculation scheme to transform 1-flat systems into that desired form. It 
should be mentioned again that we only provide sufficient conditions for a sys- 
tem to be 1-flat and that the constructive algorithm is in general not unique, 
and iterations might be necessary. Nevertheless, we believe that the presented 
normal-form is of interest in the analysis of the flatness problem, and our ex- 
amples show that this implicit triangular form can be achieved by successive 
coordinate transformations in a rather straightforward manner. Additionally 
the well known Brunovsky form for systems that are linearizable by static 
feedback is naturally included in our approach, based on proposition 14. 
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